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ABSTRACT 


Previous research of systems diagnosis algorithms have 
assumed that replacement processors are fault-free. fia 
practical applications, however, faults can occur in spare 
processors. It is shown that faulty spare processors have a 
surprisingly large deleterious effect on the speed of 
diagnosis in the universal diagnosis Algorithm_1 anausyzed by 
Smith [Ref.13]. Algorithm_1 is described as follows; 
“Replace a processor if it fails at least one test." The 
speed of diagnosis is nearly independent of the distribution 
of fault processors. That is, as long as the total number 
of tault processors is constant, the probability of repair 
is relatively unaffected by whether more faulty processors 
are in the spares or in the system. 

In this thesis we derive an asymptotic approximation to 
the probability of repair when faulty spares are present. 
An exact value can be obtained from previously known 
results. However, the calculations are extremely time 
consuming with a time complexity of order 0(4"), where n is 
the number of processors. Our asymptotic approximations 
yleld good estimates that can be calculated quickly. The 
analysis was performed by formulating the probability of 
repair calculations as a multiplication of matrices and by 


deriving approximations to the largest eigenvalues of these 


matrices. Also, faster calculations were achieved by an 


aggregation operation on the states of the system. 
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A. GENERAL 

In the past decade the concept of fault tolerance in 
digital systems has received considerable attention. This 
interest has been motivated by the need for computers that 
operate properly under hardware failures and software 
errors. Reliability is a measure of the probability that a 
system will operate without failure. Several methods have 
been used in order to achieve reliability. One method is 
through redundancy where identical units perform identical 
computations and the results are compared. Redundancy is 
quite expensive and hence is suitable primarily for the 
Situation when repair is difficult or even impossible, i.e., 
space missions. Another method is to use fault diagnosis 
and repair. In this approach the fault must first be 
detected and then the faulty component must be identified. 
However, for real time critical applications in which fast 
fault-free operation is necessary, the repair must be very 
fast and possibly not involve human intervention, i.e., 
automatic self-repair. In order to achieve automatic 
repair, either we can use spare components or we can design 
into the system a capability of reconfiguration in order to 
work with reduced capacity in the presence of the faults. 


self-repairing implies self-diagnosability, i.e., the 
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capability of automatic detection and diagnosis of faults to 
a number of replaceable nodes. Detection and isolation of 
faults involves applying tests where a test consists of 
applying a set of inputs to the system, observing the 
resulting outputs, and comparing the outputs with their 
anticipated values. 

The development of VLSI technology makes it possible to 
partition a system into replaceable nodes and the advent of 
low cost+ nicroprocessors makes possible networks of hundreds 
(or more) of interconnected nodes. The demand for high 
availability and self-diagnosability is becoming a feature 
of major importance in digital systems. Each node can be 
tested by some combination of other nodes. There have been 
several attempts to formulate repair strategies for digital 


system fault diagnosis. 


De INTRODUCING THE PROBLEM 

Previously studies of systems diagnosis algorithms have 
assumed that replacement processors are fault-free. 
However, in practice faults can exist in spare processors. 
Lee and Butler [Ref. 19] have shown that faulty spare 
processors have significant negative effect on the speed of 
diagnosis in the Algorithm_1 analyzed in Smith [{Ref. 13]. 
In this work, we develop a procedure for formulizing an 
optimal design with respect to speed of diagnosis ina 


system which consists of replaceable nodes, and we consider 
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the number of steps for repairing a system. Our analysis 
also shows asymptotic approximations to the probability of 
repair which require eigenvalues produced by the matrix 
manipulation program MATLAB. However, there is a limit on 
the dimension of matrices which MATLAB can handle. An n 
node system requires a matrix of dimension 22 x 2n, which is 


very large for even moderate size n, as shown in Table 1. 


fap Lia 


‘4 OF NODES IN THE SYSTEM VS 
TRANSITION MATRIX DIMENSION 


+ OF NODES TRANSITION MATRIX DIMENSION 

D X 
5 8 x 8 
4 16 x 16 
5 Box oo 
6 64 x 64 

10 Oo Auaeunno a 

15 B27 Ge <se7 6S 

20 1,048,576 x 1,048,576 
n billions x billions 
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C. CONTRIBUTION AND ORGANIZATION OF THE THESIS 

Many diagnostic models measure diagnosability by the 
number of faults a system can tolerate in the worst case. 
Additionally, a very conservative philosophy of system 
repair is assumed, i.e., only faulty nodes can be replaced. 
An alternative repair strategy [Ref. 13] involves inexact 
diagnosis and replacement of faulty nodes plus some nodes 
which may not be faulty. 

The goal of t**3 thesis is to derive a probability of 
repair calculation that can be performed in reasonable time. 
Previously, exact calculations were available, but required 
extensive computer time and memory. Our focus is on 
Algorithm 1 of Smith [Ret] ore In Chapter II, we survey 
the literature on system diagnosability. In Chapter III, we 
present a sequential diagnosis of single loop systems when 
Algorithm_1is applied. In Chapter IV, we derive the 
transition matrix of the probability of repair using 
Algorithm_1. In Chapter V, we present the formulation of 
our first approximation. In Chapter VI, we introduce the 
aggregation of the transition matrix PP. Conclusions and 


some suggestions for future work are presented in Chapter 


Vie 
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Ii. BACKGROUND AND NOTATION 


A. SURVEY OF THE LITERATURE ON SYSTEM DIAGNOSABILITY 

In the systems diagnosis model of reliable 
multiprocessing systems, faulty processors are identified 
from test results produced by other processors in the system 
[Ref. 1]. The goal of the diagnosis algorithm is to replace 
possibly faulty processors with spares so that the system 
consists entirely of fault-free processors. 


In all previous systems diagnosis models, it is assumed 


that, when a processor is replaced, its replacement is 
fault-free. While this assumption simplifies the analysis, 
it is not realistic. In a practical system, the most 


reliable components will be used to do the computation, 
while the least reliable ones are kept as spare processors. 
We consider the effect faulty spares have on the speed of 
diagnosis. When there are even few faulty spares, a 
Significant degragadation in speed of diagnosis can occur. 

A SYSTEM is a directed graph in which nodes represent 
processors and arcs represents tests between processors. 
Let {Uo, Urs ...Un-1 } be the set of n processors. 
Associated with each node is a status, faulty (bad) and 
fault-free (good). If there is an arc from Ui to Uj, then 
Ui tests Uj. Associated with each arc is a test outcome, 


which is generated as follows: 


ae 


The outcomes of the tests are represented by binary 
values on the arcs (aij) where aij is ’O’ if Ui evaluates JU; 
to be fault-free, and is ’1’ if Ui evaluates Ujto be faulty. 
If Ui itself is faulty then the evaluation of Uj is 
unreliable, Since aij can assume ’O’ or ’1’ regardless of 
the status of Uj. If Ui is fault-free, its test results are 
correct; it passes fault-free processors it tests and fails 


faulty processors. 


aij = 1 if ™ evaluates Uj as faulty, 
and 
ais ih =a if Ui evaluates Uj as fault-free. 
Figure 2.1 summarizes these assumptions. Here a open circle 


represents a fault-free processor and an asterisk represents 
a faulty processor. For example, the topmost arrow between 


open circles represents a test by a fault-free processor Ui 


of a fault-free processor Uj. It follows from the previous 
discussion that the test outcome aij is 0. This is shown in 
Figure 2.1. 


Figure 2.2 shows an example of a five processor 
multiprocessing system where each processor is tested by 
exactly one other one. It is called a single loop system. 
In Figure 2.2 the complete set of test outcomes is shown. 
Such set is called a SYNDROME. 

The object of a diagnosis algorithm is to determine 
Which nodes are faulty given the syndrome. Because of the 


arbitrary test results by faulty nodes, it may be impossible 
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ao fault-free node 
* faulty node 


Figure 2.1 Assumed test outcomes in the 


Preparata-Metze-Chien model 





Figure 2.2 A five_processors single loop 
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system and 


to uniquely identify the faulty nodes. A step in the 
algorithm consists of idenmiaeeine nodes to be replaced, 
replacing them with spare nodes, and generating the next 
syndrome. The diagnosis is complete when all test results 
are 0. It is assumed that no fault-free node becomes faulty 
during diagnosis. 

Figure 2.2 shows how the presence of faulty processors 
inhibits the identification of faulty processors. For 
example, in Figure 2.2 there are thre- fault-free processors 
(open circles) and two faulty processors (asterisks). If 
both faulty processors produce pass (0) test results we have 
a syndrome that is also produced by the set shown except 
that U2 is fault-free. It is not known beforehand exactly 
which processors are faulty, and there is an ambiguity. We 
could resolve this by replacing U1 only and applying tests 
again. Since after this U1 is fault-free, we can identify 
U2 as faulty and replace it, thus completing the diagnosis. 
However, this requires two diagnosis steps. 

Mallela and Masson [Ref. 9] have considered a system in 
which faults are intermittent. That is, at the time of one 
test application a tested node U may be faulty, while ata 
later time it may be fault-free. Butler [Ref. 16] developed 
the relationships between diagnosability of general systems. 
That is, three models are compared; one model is the 
conventional system with binary-valued test outcomes, the 


other has three valued test outcomes, (where the third value 
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is missing test result) and the third model corresponds to 
permanently and Seoeettt ent Ie faulty processors. Also 
Butler [Ref. 15] extended the technique to the case where 
intermittently, as well as permanently, faulty processors 
are presented. Karunanithi and Friedman [Ref. 11] 
considered replacement algorithms for the Preparata-Metze- 
Chien model in which fault-free nodes could be replaced. 
All previous studies assumed that no fault-free node could 
be replaced. Chwa and Hatimi [Ref. 14] analyzed in detail 
diagnosis under the conditions proposed by Karunanithi and 
Friedman [Ref. 12]. In both studies, the goal of the 
diagnosis strategy is to determine the smallest set of nodes 
which contain all faulty nodes and perhaps some good nodes 
as well. Maheshwari and Hakimi [{Ref. 9] considered a 


probabilistic approach for fault diagnosis. 
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TTI. SEQUENTIAL DIAGNOSIS 


A. k_ STEP t/s_FAULT DIAGNOSABLE SYSTEMS 
A new measure of system diagnosis, t/s diagnosability 
originally proposed by Friedman [{Ref. 6], is used to study 
the diagnosability of digital systems. Friedman used this 
measure to study a canonical class of systems, “Single Loop 
Systems.” Generally, the system is repr--ented and analyzed 
by means of the graph theoretic model of Preparata et al. 
fRef. 1]. 
Let S be a single loop system, a system which has 
exactly one test link outgoing from node Ui, i =0, 1, 2 
., n-1l, which is connected to node Uit+l( mod n) where nis 


the number of nodes, as in Figure 3.1 





Figure 3.1 A single loop system 
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The: question we wish to resolve is how this system 
should be repaired. We assume that the probability that i 
nodes are faulty is less than the probability that it+1l nodes 
are faulty. 

Friedman defined a new measure of system diagnosability 
related to this concept. 

Definition 1: A system S$ is k_step t/s (t-out-of-s) 
diagnosable if and only if by no more than k applications of 
the diagnostic tests any set of nf < t faulty nedes in S can 
be diagnosed and repaired by replacing at most s nodes. 

Sayrously, so > t and n > ‘s- If s = n, repair is 
trivial since the entire system is replaced. A measure of 
the diagnosability of the system is the average or the 
maximum value of s-nf over all system fault conditions with 
nf < t faulty nodes. 

Here we will introduce three replacement algorithms 
which have been presented in Refs. 1, 12, 13, and 18. 
Algorithm_1 and 2 are similar in Smith [Ref. 13]. The 
algorithms are as follows: 

Algorithm. 1: At each repair step, perform the tests 
and replace the nodes which fail at 
least one test with randomly selected 
nodes from spares. Replaced nodes are 
placed back into the set of spares. If 
all test results are Q (pass), then 


the system is assumed to be correct. 


a3 


Algorithm_2 : At each repair step, perform the tests 
and replace the nodes which fail the 
maximum number of tests. Replaced 
nodes are placed back into the set of 
spares. If all test results are QO 
(pass), then the system is assumed to 
be correct. 

Algorithm_3 : At each repair step, perform the tests 
and replace the nodes which fail +*+4 


—— 


maximum number of tests. i¢ the 
number of modes in Spare-1l is not 
enough, randomly select any additional 
needed spare nodes from opare-2. 
Replaced nodes are placed back into 
the Spare-2. [If all test results are 
O (pass), then the system is assumed 
to be correct. Initially, all spare 
units are in Spare-1, and Spare-2 is 
empty. 

A good sequential diagnosis strategy will identify many 
faulty nodes for replacement at each diagnosis step. It has 
been shown (Refs. 13, Seed) that Algorithm_1 is 
Significantly faster than Algorithm _2 when spares are fault- 
free, while the reverse is true when some spares are faulty. 
Algorithm _ 3 is faster than either Algorithm_1 and 


Algorithm_2, in the presence of faulty spares; it was 
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designed to compensate for the negative effect of faulty 


spares. 


Be SEQUENTIAL DIAGNOSIS OF SINGLE LOOP SYSTEMS USING 

ALGORITHM_1 

In this work, we selected Algorithm_1 because it 
results in the most tractable analysis. We assume symmetric 
invalidation [Ref. 1] where the test outcome of a test by a 
faulty node is 0 or 1. The possibility of exch occurrence 
1s assumed to be 1/2. 

Repair means the replacement of possibly faulty (bad) 
nodes by a set of spare processors. For sequential 
diagnosis, every time that a repair step has been terminated 
all tests will be performed to check whether the system is 
correct or not. It is assumed that, for more practical use, 
even if both nodes Ui and Uj are left untouched during a 
replacement process, Ui will be tested again at the next 
step because some permanent faulty nodes could become 
meaty . 

As mentioned before, Algorithm_1 is: "Replace a node if 
it fails at least one test." <A repair step in the algorithm 
consists of identifying nodes to be replaced, replacing them 
by spares some of which are faulty. All replaced nodes are 
placed back into the set of spares since they might contain 
fault-free nodes. We do not assume that the new replacing 


node is fault-free. That is, there can initially exist 


ZS 


faulty units in the set of spares. This situation is shown 


VY Eee see oe An arrow trom the spares to the system 
represent a choice of spares to use as replacement 
processors. An arrow from the system to the spares 


represents the placement of replaced processors back into 
the set of spares. 

Another problem created by the presence of faulty 
spares is a degradation in speed and efficiency of 
diagnosis. In the worst case, if only faulty units “re 
involved in the replacement, then an infinite cycling of the 
procedure is possible. The problem can be solved by having 


sufficiently many fault-free spares. 
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Figure 3.2 Sequential Diagnosis of single loop 


systems using Algorithm_1. 
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IV. CALCULATION OF THE PROBABILITY OF REPAIR 


A. METHOD OF APPROACH 

We seek the probability of producing a completely 
fault-free system after the application of the Algorithm_i 
given that there is a fixed number of initially faulty nodes 
in the system and a fixed number of faulty spares. Because 
of faulty spares, on each application of the Algorithm_1, 
faulty nodes can be introduced into the system. Further, 
Algorithm_1 may incorrectly assess the status of a node and 
so fault-free nodes may be replaced by faulty nodes. 
Therefore, it is possible that an application of the 
algorithm can produce a net increase in the number of faulty 
nodes. In calculating the probability of repair [Ref. 18], 
we assume: 

1. Any arrangement of a fixed number of initially faulty 
nodes is as likely as any other arrangement. 

2. Each test by a faulty node is pass or fail with 
probability 50%. Test results are generated at each step of 
the Algorithm_1, so that the test result by faulty node Ui 
of node Uj may not be the same after the application of the 
Algorithm_1 as before, even though neither Ui nor Uj is 
replaced. Such behavior is a characteristic of intermittent 


hPaults. 
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3. Spares are chosen randomly, and the total number of 
spare processors equals or exceeds the number of system 
nodes identified for replacement. 

Let Fi be a set of faulty nodes in the system after the 
i-th application of an Algorithm_1, with Fo denoting a set 
of faulty nodes just prior to diagnosis. Let pi(F —> F’) 
be the probability of attaining fault set F’ afteri 
applications of Aizen to F. We seek the probability 
of repair, P (nf, i), the probability that, starting ag 
fault sets of size nf, there will be no faulty nodes after i 


applications of Algorithm_1. That is, 


P(nf 1) = > P; (F N) (eqn. 4.1) 
| L | 
ce 
' "0 — No: emoty (null) set 


where the sum is over all fault sets of size nf and where N 


1s the empty (null) set. 








Pinfia) => > D . F_ ptr. _ +N} 
iis i-] 
| Fy =nf Fy 
(eqn. 4.2) 


where the left sum is over all initial sets of nf faulty 
nodes and the right sum is over all sets of faulty nodes. 
The derivation which follows is for p (Fi-1 —> Fi), where 


Fi-1 and Fi are any faulty § set. The second factor in 


aS 





Equation 2 is then obtained directly by a substitution of N 


for Fi, while the first factor is obtained iteratively. 


We derive p (PFi-1 — > Fi) by tracking the number of 
nodes which undergo various transitions. We consider the 
derivation of p (Fi-1 —> Fi) for a two node system using 


Algorithm_1 in next section. For all cases, the approach is 


to derive two probabilities. That is, 
ej) = 2 PLE 7 alpla—+F) 
a 
(eqn. 4.3) 
where p (Fi-1 —— > a)is the probability Fi-1 produces 
syndrome a, p (a ——> Fi) is the probability a choice of 


Spare nodes, as determined by "a" and Algorithm_1, will 


produce Fi, and the sum is over all possible syndromes “a”. 


We have, 

na fC nt 

o(F ——»a)=2 (eqn. 4.4) 
i- | 

where nf = | Fo], the number of initially faulty nodes and 

nf’ = (ine9 |. ie miuMivemEOlm fawmey nodes alter the “i-l th 





application of the Algorithm_1. p (a > Fi) is just the 
number of ways to choose spares s0 that Fi is produced 
divided by the total number of ways to choose sparés. 

While we show derivations for certain 2, 3, 4, and 5 


node single loop systems, the derivaticns apply to all 


designs. 


Be A MATRIX REPRESENTATION OF THE PROBABILITY OF REPAIR 

CALCULATION 

We show in this section that the probability of repair 
as calculated in Section A, can also be computed using 
matrix operations. This not only gives additional insight 
into the problem, but also facilitates a calculation of 
approximate values Of the probability of repair. 
Approximate values are necessary when the number of nodes is 
such that the calculations (either by Equation 4.3 or Mave. 
manipulations) require inordinately large amounts of 
computer time and memory. For example, the calculations by 
Lee [Ref. 19] requires computation times on the order of 
days for even moderate sizes of nodes (n = 15). Approximate 
values are calculated in Chapter JV. In this section, we 
illustrate the matrix implementation of Equation 4.3 by the 
two node system. 

i One application of Algorithm 1 

When the test outcome 15 252e70 (aba = Bape 

Algorithm_1 replaces one good node "b" from 8 good spares. 
The contribution of this test result to Pi (Fo—>N) is (1/2) 
(8/10), since there is a probability of 1/2 that aba = 0 and 
there are 8 out of 10 ways to choose a fault-free spare. 
This is shown on the next page. When test outcomes one (aba 
= 1), Algorithm_1 replaces two good nodes "a" and "b" from 8 
good spares. The contribution of this test result is (272 


28/A, where A is the total number of ways to choose spares. 
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Here, 1/2 is the probability aba = 1, and 28 is the number 
of ways 2 good nodes can be chosen from 8. A is calculated 
as shown on these page. Here randomly chosen spares could be 
all good spares with one possible way or all bad spares with 
one possible way or one good, one bad spare with two 
possible ways (A= 61). We calculate the probability of 


repair after the first application of Algorithm_1 as 0.6294. 
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Figure 4.1 Illustration of First Application (i = 1) 


Using Algorithm_1 for Two Nodes System. 
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To do this we 


need to know the number of ways to choose k 
nodes from n. This is, 
J 
: + nm! 
F a ee (eqn. 4.7) 
K 


(n-k)Ik! 


Let ipj be the probability of transition from state i to 
state j, 


where a state 18S a specification of which nodes are 
“aulty and which are fault-free. Figure 4.1 shows the 
possible transitions between states, 


assuming initially that 
nNodewis Pauli 
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Figure 4.2 Illustration of Secend wapelicae1en i) ee 


using Algorithm_1 for two nodes system. 
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a. Two a mea tems Or A Ler. thm 
Now let’s calculate the probability of repair 
after two applications of the Algorithm_1 for a two node 


system. 
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We calculate 1lpO = 0.6294 which is the prebability of remem 
from state 1 to state 0. Also Op0O = 1.0 since Algorithm_1 
leaves the system unchanged if it is in state QO. 
The calculation of ipl: 

When test outcomes aba = O, Algorithm_1 replaces one bad 
node "b" from 2 bad spares. When test outcomes aba = 1, 
Algorithm_1 replaces one good node "a" from 8 good spares, 
one bad node “b" from 2 bad spares. As shown in Figure 4.4, 
the result is lpl =s@mZgiz- 

S The calculaticnees. 1pZ: 
When test outcomes aba = 0, the probability of transition is 
O because Algorithm_1 does not replace node “a”. So node 
“a’ does not change from good to bad. When the test 
outcomes aba = 1 Algorithm_1 replaces one good node “b" from 
8 good spares one bad node “a from 2 bad spares. The 


result. 21s. 1 p2 (= 30. 


he calewiatiaan of 1 ps: 
When test outcomes aba = 0, the probability of transition 1s 
O (likewise calculation of 1p2). When test outcomes aba = 1 
Algorithm_1 replaces two bad nodes “a’ and "“b'" from 2 bad 
spares. The result is 1p3 = 0.0082. 

The calculation of 1p0: 
When test outcomes aba = O Algorithm_1 replaces one good 
node "b" from 8 good spares. When test outcomes aba = 1 
Algcrithm_1 replaces two good node "a" and "b" from 8 good 
spares. The result is 1p0 = 0.6294. 

The calculation of 2p0: 
The result is same as 1p0O. Since state 1 and state @2 are 
the same except for a rearrangement of faulty nodes. 
(Assumption 1). Therefore, 2p0 = 1p0 = 0.6294. 

The calculation of 3p0: 
Unlike the cases discussed previously, Algorithm_1 may fail 
to replace any node. That is, when both faulty nodes 
produce pass test outcomes Algorithm_1 makes no replacement. 
When aba = aab = 1, the transition to state 0 requires that 
Algorithm_1 replace both units with fault-free spares. This 
is with probability 0.1666. When aab = O and aba = 1 or 
vice versa only one faulty node is replaced. For these test 
outcomes attaining state 0 is impossible. Thus, the total 
probability of transition from state 3 to 0 is 0.1666. The 


number of ways, B, to make choices from spares is 54. 


35 


os Three applications of Algorithm 1 

Now let’s consider three applications (i = 3) of 
Algorithm_1 for two nodes system. 

We can realize that this calculation will be just 
a repetition of the calculation in Figure 4.2. This is 
exactly same transition from one state to another state. 

We already calculated first row and second column 
of transition matrix P before. Note Opl = Op2 = Op3 aaa 
because if the *~itial state is 0, the final state will also 
be O. Because of symmetry lipl = 2p2, 1ip2 = 2pl1, 1ip3 = 2p3 
and “Sieote=—30c- (Assumtion 1) 

The calculation of 3pl: 

Here nf’ = 2 when test outcomes aba = O. There is no 
possibility, because Algorithm_1 does not replace node ‘a" 
from bad to good. When aab = O and aba = i, Algorithm_i 
replaces one good node "a" from 9 good spares. When aab = 
aba = 1, Algorithm_1 replaces one good node "a" from 9 good 
spares, one node "b" from 1 bad spares. The result is 3pi = 
0.2667. 

The calculation of 3p2: 

The result will be same as 3pi1. This is another arrangement 
Gt S3p ie There is no difference (Assumption 1). So 3p2 = 


Spl =" Onze6 (o 
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Figure 4.3 Illustration of third application (i = 3) 


using Algorithm_1 for two nodes system. 


The calculatrzonsot. sco. 


When test outcomes aab = aba = 0, there is no replacement 
and the probability is 1.0. When aba = 1 Algorithm_1 
replaces one bad node "“b”" from one bad spare when aab = 1 


Algorithm_1 replaces one bad node “a” from one bad spare. 
When aab = aba = 1, the probability of transition is 0. 
Because the system has only one bad spares. Therefore 
Algorithm_1 does not replace two bad nodes “a”"” and "“b” from 


one bad spare. Tre result is 3p3 = 0.3. 
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NemiGmrecdmEtiat seacae  arpplication of Algorithm_l 
has the same transition states. Therefore the P matrix is 


Seeds a transition matrix of the probability of repair. 
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V. CHARACTERIZATION OF APPROXIMATION USING MATRICES 


A. FIRST APPROXIMATION 

While the matrix multiplications for the examples in 
Chapter IV can be done easily by hand or computer, the 
multiplications increase dramatically as the number of nodes 
increases. For an n node system, the transition matrix 
dimension is 2" x 2M as shown ir Table 1. Its size is very 
large for even moderate values of n. 

We illustrate our method of approximation using the two 
node system. For example, if we have two nodes (n = 2), we 


get four states (22) indexed by s = 0, 1, 2, 3. 
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Figure 5.1 Four states to repair the two nodes system 


Let Ps be the probability of being in state s before 
diagnosis and Ps’ is the probability of being in state s 
after diagnosis. Therefore, after one application of 


Algorithm_1, we have, 


DT Se (eqn. 5.1) 
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Now we seek, 


px = Xx 


Boe S ep EOSabl  ityeot transition 


applications of the diagnostic algorithm. 
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Markov chains provide one example in which we are 
interested in computing PKX. The transition matrix P has 
the property that all entries are nonnegative and the sum of 
the entries in any column is 1. X is the probability 
transition vector before diagnosis. Then PX is the 
corresponding probability transition vector after one 
application of the diagnostic algorithm. Similarly, P2X is 
the probability transition vector after two applications of 
the diagnostic algorithm and, in general, PKX gives the 
probability transition vector after k applications of the 
diagnostic algorithm. 

In order to derive a simple closed form approximation 
to the probability of repair, it is necessary to derive the 
eigenvalues of the transition matrix. 

Note: Suppose that P is an nxn matrix and V is a column 


vector with n components such that; 


Oo Y=E ¥ (eqn. 5.5) 


for some scalar E, where E is eigenvalue, V is eigenvector. 


It is an easy exercise to show that; 


kk kk 
DW Y (eqn. 5 769 


Thus PKX is easily computed if X is equal to this vector V. 
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i. Diagonalizing a matrix 

A square matrix is called diagonal if all 
nondiagonal entries are zero. Later we will indicate the 
importance of being able to compute PKX. We show that if P 
has s distinct eigenvalues, then PK in this computation of 
PKX can be essentially replaced by Dk, where D is a diagonal 
matrix with the eigenvalues of P as diagonal entries. It 
can be seen that DK is the diagonal matrix obtained from D 
by raising each diagonal entry to the power of k. 

THEOREM: (Diagonalization) 


Let P be an nxn matrix having n distinct 


eigenvalues Ko, Fi, H2, ..., Es. Let Vo, V1, V2, ..., Vs be 
column vectors in Ra such that Vi iS an eigenvector 
corresponding to Ei for i=0, 1, 2, ..., s. Set C be the 
nxn matrix having V as i-th column vector. Then C is 


invertible and C-1 P C is equal to the diagonal matrix D. 


D = | eed (eqn. 5.7) 

— , 
f 
|} 
LY E 

(ame, 2 a] 

EeOor : 
We know that, 
? ¥, = EY. (eqn. 5.8) 
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[0 PC = CD and C7! FP € = 2 @We ss 


Po | ot Ty ffoe 
i : | O) 
P| 

/ 

: 





Now we snow how the computation of PK can be essentially 
reduced to the computation of Dk. Let P , Ei, Vi, and C be 
as in Theorem. Then PK = C Dk C-!1 for each positive integer 


k. From C-! P C = D, we Gbtain P = Cee. heuer 


K ae =: ae 
P'=(CDC “CDC -)......... (GDE =) (eqn. 5.10) 
—$_$§_——— _ k factors ————e 
Clearly the adjacent terms C-! C eance We get 
Pk = C DK C-1. Then we say that P has been diagonalized by 


O 
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An nxn matrix with n distinct eigenvalues is 
) Similar to a diagonal matrix. It is not always essential 
that the eigenvalues be distinct. 


ae, Computing PkKX 


Let P be a diagonalizable nxn matrix and let Eo, 


) re, SCti«<i2,, ; Es be the (not necessarily distinct) 
| 

eigenvalues of P. het Vo, Vi, Y2, ..., Vs be a basis for 
) Rn, where Vi is an eigenvector for Fi. We have seen that if 


C is the matrix having Vi as the i th Some oti at en 


‘a e 
. a } | 
\ 
-1 : | i | 
ose = ff = a (eqn. 5.11) 
i ’ 
| 
| LY E | 
es S 
tie 


Let X be any vector in |!R".. We show that from Equation 5.6 


we get, 
[_ > 
| | | i 
! i | 
| | i 
K k -] K ke 
4 A = C ry = | ‘a Cc ee - K 8 ma 
A C x Ey On a W str Gone w carole ra a tate iw lems’ s S| C x 
| : | 
| | 
eS 
(eqn, 5,12) 
Now C-1 X is a column vector. It may be any column vector 


ar Rn. We set, 
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C X= Cilia OC (eqn. 5.13) 


c 
o 
fi 
[= 


Then, we form 
oky og cky ‘+g eX + + Gd EY 
- “9"o i fie $$ $_ (eqn. 5.14) 


This equation expresses PkKX as a linear combination of the 
eigenvectors Vi. X is a vector whose components are the 
probability of being in a particular state of the system.- 
The probability vector at the next stage of the process is 
found by multiplying the present probability vector on the 
left by a transition matrix P. Tllustrations of Gigs 
Situation are provided by Markov chains and the generation 
of the Fibonacci sequence. 


We are interested in the long term outcome of the 


process. That is, we wish to study FKX for large values of 
k. In particular, suppose that (|Eo| > |£1 | ) so that E 
is the unique eigenvalue of maximum magnitude. If do f 0 
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and k is large, the vector PKX is approximately di: Ei1* Vi in 
the sense that || PKX - do Eok Voll is small compared with 
Pk X II 

Let’s consider our example again. We can get 
@igenvalues and eigenvectors from MATLAB computer program 
for two node system. PC-MATLAB is a general program for the 
scientific and engineering numeric calculations. The name 
MATLAB stands for "“MATrix LABoratory"”. PC-MATLAB requires 
the following hardware and software: 

Sie Pee FC/XL, FPUC/ATS or compatible MS-DOS 

computer. 

- At least 320K of memory. 

mows version 2.0 or higher. 

=weeo numeric coprocessor chip. 


- At least one DSDD floppy disk drive. 





; -J.00090 1.0000 1. OOOU 
Gee +1 0000 -O4635 -Oo215 (eqn. 5.15) 
| WUC s-O 66s 0.52 15 
L Oi ee an 
Ee 
E soa ee | 
. 5 
a Se | 
| 0.1000 | (eqn. 5.16) 
D - 0.4043 | 
7 oy ny 7c 5 ; 
jee eer 
(a. aw) 
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| 1.0000 1.0000 = 1.0000 
ct_] 0 +0.5000 -0.5000  -9.0000 


(eqn. Saaz 
O -0.76984 -0.7694 -3,.9225 
| =U. 250, =O.2ee 2.2335 44 


Substituting these eigenvalues into Equation 5.14, 


K 


K K K K 
P -Aa=d (i ¥e+d (0.1 Y +d (0.40435) V d_ (0.25d1 Y 
gi ~G OG a 3) 7 °dt0:2 ) ¥. 
(eqn. 5.18) 
Then, 
K 
P X-d _¥ 


K K K 
= ’ + f 2 + 7 
oo S (0.1) we d (9.4043) Y, d, (0.25538 Y 


(eqn. 5.19) 


As k increase, the middle term on the right dominates and as 


an approximation, we have 


"tO 
> 

| 
cL 


(eqn. 5.20) 
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From our definition 


foci cs) 





where X’ vector is the probability of being in each states s 


after diagnosis. 


a | 
aA = d (eqn. 5.22) 
— “ and = 
f 4.0000 10000 1.5000 P J 
| 0 9.5000 -9.5000 -0.0000 [1 Fy qd 
| = | (Ceadnn Sac3) 
O -O.7394 -0.7334 -3.3556 | PA de 
i? | a 
| 9 “G.2590 <O.2505 2.59555 | Ps | 9 
We have, 
d.=P.+P +p. + ~ 
Q O | 2 PS =| (eqn. 5.24) 
Also, 


ean. 5.25) 








and arranging our Equation 5.10. We get, 


ee 
| 0 

P 0 k 

a = d,(0.4043) ¥ (eqn. 5.26) 

FY 0 > 2 

> 

' lo 
| ~~] Le 


We get d2 from Equation 5.23 and V2 third column vector from 


Crmatrix. 256, taace 


d. = 0.0 Fo ~0.7694P, ~0.7694P., —5.9555R, (eqn. (5am) 


Hoon 
= | “0.46726 (eqn. 5.28) 
2 eae 
| -0 3729 | dé | 
We define, 
Pp =i PP =P =P = ——® tst column cf p* 
0 1 2 S a 
P =i = = a eS =] nmr a 
i Pa e Po OO ——2} End column oi P 
D =1 Pp =P _PD =O ——& Grd column of p* 
2 0 oA = ee an 
P. =1 Pa ee 17 P, =O ——@® {th colurmn of P 
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Then we find all 


=| -0.7694P -0.7694P -s933er, 
L | 2 


«hx | 











— as 4 mm = K 
UG50.(0 s0d5, 

= = 7 ¢ ..X 
0.555710 A085) 


0.055 110.2045) : 


computations of Pk 


v k 
“7. FO9dU dS) 0. 7SE9KO dOdT) - F95T865 was) 
_ = 7B sm _.A 
0.555 .(G.—0a5) 
toler K 
O.SE57(0 WaZ) 


0.05501 (0 40a)" 


1.G060 
kK =) d636 

(0 4047) 
Ixa -0.5656 
=). 0729 


<4 


(eqn. 5.29) 
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x 
we ee 


* ts. 
_ = £ am e = 
coy ¢ (O..O42) 


k 
|) 828710. 4042) 
ee) 
mo 


(eqn. 5.30) 


a 





on ae i ke k 
| ae f ie Tess O jp le PSSM E ) l= Pewee) | 
as | a 
: ten K ere 
Poy | 0 OrS67(E) O5S670E) Ween. )e 
K I 1 1 | 
> = 
: ek re eee ime i 
pP 0 OSSS.1E ) 0 so-0e,) Peete.) 
2 | | 
eae ; K eee re. 8 
| Ps | | ° OUSSiKE )  OUSSI(E 0.2560 2,1") 
ws 
a (eqn. 5.31) 
k This vector ts orobsoilitu of transition states arter kK 
soolications of the diagnostic sidgomtnm |. 
e ° a @ 8 _ = K = = = %S 
P 2: 1.0 (P -# =e ae H-0 $69 cE iy iP )- 7.925 55¢ eae > 
Q Det ae } ie 1 oo 
P* =P ° SISSTE UP O ye! 3 ‘o 
= =a RP uy = = ° Ss 3 e a 
' 4 ov Vos , 2 ( ye 
Pee OUso lc sip +7 )*0 S5G0l & b 3 
See a) ee (eqn. 9.32) 
=N aNd7 Second largest 2igenvalue 


E 
| 


on 


Be COMPUTING SECOND LARGEST EIGENVALUE BY THE POWER METHOD 


As shown in Equation SO. the second largest 
eigenvalue must be calculated. Several algorithms have been 
developed for calculating the largest (or smallest) 


eigenvalue, and no one method can be considered the best for 
all cases. Computation of the eigenvalues of a matrix is 
one of the toughest jobs in linear algebra. We use the 
program MATLAB in this thesis. Tere are other methods; for 
example, the power method, Jacobi’s method, and the QR 
method. 

We also use the power method. It is especially useful 
if one wants only the eigenvalue of largest (or of smallest) 
magnitude, as in many vibration problems. So we need second 
largest eigenvalue here. We know that our first largest 
eigenvalue is always 1. 

Definition: An eigenvalue of a matrix P is called the 
dominant eigenvalue of P if its absolute 
value is larger than the absolute values 
of the remaining eigenvalues. An eigen- 
vector corresponding to the dominant 
C€igenvalue is called a dominant eigen- 
VECLoOrRRer Cs. 

Since P is diagonalizable, there exists a basis {Vo, 

Vig). .2e -¥saanor !Rn composed of eigenvectors of P. We 


assume that V1 is the eigenvector corresponding to E1, and 


D2 





mide the numbering is such that Ko.) ey E: | 2 E2 | 
2 Pees > | Es | . Let W be any nonzero vector in | Rn. 
Also the sum of the values in W is 1. Then we can 
approximate an eigenvector of P corresponding to the 
dominant eigenvalue £1 by multiplying an appropriate initial 
approximation vector W repeatedly by P. 

Repeated multiplication of W by P may produce very 
large (or very small) numbers. It is customary to scale 
after each multiplication to keep the components of the 
vectors at a reasonable size. After the first 
multiplication, we find the maximum m of the magnitudes of 
all the components of PW and apply P next time to the vector 
W1 = (1/m)PW. Similarly, we let We = (1/m )PW1, where mi is 
the maximum of the magnitudes of components of PW1, and so 
on. We summarize the steps of the Power method in the 
following Table 2. 

TABLE 2 


THE POWER METHOD FOR FINDING THE SECOND LARGEST 
EIGENVALUE OF P 


Step 1 Choose an appropriate vector W in 
(Rn as first approximation to an 
eigenvector for fi. 

Step 2 Compute PW and the Rayleigh 
quotient (PW . W)/(W . W). 

Step 3 Let Wi = (i/m)PW, Where m is’ the 
maximum of the magnitudes of 
components of PW. 

Step 4 Go to step 2, and repeat with all 
subscripts increased by 1. The 
Rayleigh quotients should 
approach El, and the Wj should 
approach an éigenvector of P 
corresponding to £1. 
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We find second largest eigenvalue for 2, 3, 4 and 5 
nodes system by power method in the following pages. 

Note that if P has an eigenvalue of nonzero magnitude 
smaller than that of any other eigenvalue, then the power 
method can be used with P-i to Girma this smallest 
eigenvalue. Recall that the eigenvalues of P-i are the 
reciprocals of the eigenvalues of P, and the eigenvectors 


are the same. 


4 


selected arbitrary matrix W to get second largest 


@igenvalue of the (1/2 nodes bad and 2/10 spares bad) 


system. 
1 
= 1 
ARBITRARY VECTOR W = (eqn. 5.33) 
i 
=a 





TABLE 3 





POWER METHOD FOR TWO NODE SYSTEM 


Rayleigh Quotient 


1. step of power method OF 3334 

ao / O22 5256 

S: . 0.4718 

4. °" 0.4423 

Ss «OF 0.4267 

oC 0.4179 

(ae OrairZe 

8. 0.4096 

9. " ‘i 0.4077 

nO. ii iu 0.4064 

im bs 0.4057 

12. ‘i is Om4952 

i... 0.4049 

14. “ 0.4047 

i. °° 0.4045 

| we." i 0.4044 
it, i ie 0.4043 ->second 


largest eigenvalue 


After 17 steps of the power method, we get our second 


largest eigenvalue. 
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Selected arbitrary matrix W to get second largest 
eigenvalue of the (2/3 nodes bad and 2/10 spares bad) 


system. 


al 
ARBITRARY VECTOR W = fl (eqn. 5.34) 


8x1 


TABLE 4 
POWER METHOD FOR THREE NODE SYSTEM 
Rayleigh Quotient 


step of power method 0.2309 
ss : 0.5084 
S664 
poclsiil 
OalS 
1 206 
251 
one): 
ovo O 
ES 
p5ao 1 
.5150 -> second 


moon Ddonh WHF 


0 
Q 
0 
Q 
0 
0 
0 
0 
0 
0 


largest eigenvalue 
Note: After 12 steps of the power method, we get our second 


largest eigenvalue. 
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Selected arbitrary matrix W to get second largest 


@igenvalue of the (3/4 nodes and 2/10 spares bad) system. 


ARBITRARY VECTOR W = im (eqn. 5.35) 


TABLE 5 
POWER METHOD FOR FOUR NODE SYSTEM 
Rayleigh Quotient 


. 2004 
.7859 
. 1034 
Poo 52 
“Soe 
soca 
HO24iG 
6256 
fO204 
.6232 -> second 


step of power method 


OOOND OSH WMH Fr 
OOo 0o00000 0 


largest eigenvalue 





Note: After 10 steps of the power method, we get our second 


largest eigenvalue. 


Sy if 


Selected arbitrary matrix W to get second largest 


eigenvalue of the (4/5 nodes bad and 2/10 spares bad) 


system. 
iL 
1 
ARBITRARY VECTOR W = 1 (eqn. 5.36) 
= 1 
1 
ay 
6x1 
TABLE 6 
POWER METHOD FOR FIVE NODE SYSTEM 
Rayleigh Quotient 
1. step of power method 0.1626 
Z 2 0.8461 
3 0. 6506 
4 i 0.7986 
5 0.7605 
6 0. 7401 
7 0.7305 
8 . " i 0. 7282 
Guay 7 O07 2a3 
10 . : 0 Jf Za 
11 0. 72S 
152 0 .. F230 
13 0.7229 
14 0.7228 ->second 


largest eigenvalue 


Note: After 14 steps of the power method, we get our second 





largest eigenvalue. 
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VI. AGGREGATION OF THE TRANSITION MATRIX 


Already we have transition matrix P for two nodes 


system which is: 





| 60 IpQ 2po Spu 1 0.6294 0.5294 a 
Ool tIpt 2p! spl OMe fon Oe U.256 
: Mee Pe 282 op2 / Ome Urrre ore |= OU 2 ooe 
OOS ips 295 Eos lgng | YV OONTS O.NOG2 0.3000 laws 
(ane GL) 


We can illustrate our original transition matrix P as 


in Figure 6.1 





23p3 Ze 


Figure 6.1 Illustration transition states of P 


oo 


The asymptotic approximation to the probability of 
repair requires eigenvalues which in turn are produced by 
our matrix manipulation program MATLAB. However, there is a 
limit on the size of matrices which MATLAB can handle. With 
the probability of repair problem as stated, ann node 
system would require a matrix of dimension 2" x 21, which is 
very large for even moderate size n. To handle larger 
systems, we choose to represent transition probabilities in 
a different way which is called AGGREGATION. 

Let P be the transition matrix between specific states. 
Let Pa be the transition matrix between aggregated states. 
We seek the elements of Pa as a function of the elements of 
P. Consider Figure 6.2, which shows the transition between 


all states in the system with aq (initial) faulty nodes to 


| 
or 


ae ‘~ 


O-.Q0 





all states with p (final) faulty nodes. 
| 
| LoS. ee 


Ca 
(pI GD) 


Og 


QO 


a W@W C 
i; Sice c n 
| ees ~ cae p 
’ i 
System has System has 


q faulty nodes p faulty nodes 


Figure 6.2 Aggregation of system 
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= s n 
— — q 


eqn 5.2 ) 


TO Summation af all possible orobabilities fram initial 
———° state with a faulty nodes to final state with p faultu 
nodes. 
— summation of probabilities from one state with g faulty 
2 * nodes to another states with p faulty nodes over all 
possible probabilities. 
Let’s go back to our example for two nodes system. Now 


we calculate that aggregated transition probability. 





12n12 4 rs 
A -— f°} Fa 
rf =Z 
——$— — 
_ /_* 0 : 
0 sai ee oe 
fon } | 
Ag fo? 
‘if _——— 
* §=2 = oe 
| wee os YEO eo 
a yy >} a 
“hx , B 
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ip izip2e2pi+202 025 12-01 312"0 Ig ieee 











120 t= = =0.3624 
: 2 | a 
; (eqn. 6.2) 
1200: 
a =. 0 * 
( n Se ( n 5 
q aa; e- p 
<4 
1p942n0 9 UWoesd+u 6294 
12p0= > = SO=0.6294 (eqn. 6.3) 
“\ 
Lg - 
1233 
(7) = =e 
q <a + p 
$— 
lo3+203) 8690 0082+0.0082 
I2pZ= = = ————_———._ -0 0082 (eqn. 6.4) 
4g “ 
(eqn. 6.5) 


——S a a eee 


a =} 
i; OOO izZpo ial P 0.6294 0.1666 


! 
! | 

Pex Opi2 i2pi2 3soi2i-! q 03624 0.5334 
se 
cor ae a: va 0.6082 0.5000 |... 
2 at ea 


0012: 


ey 
fo 2) 
q | oe ee 
it. 











Op !+0p2 0.0 +0.0 
Op12= = =0.0 (eqn. 6.6) 

| | 

g 
SoD be 

0 = 

oe are = ~ (° } 

q — <n : 


Spi+3p2 0.2607+0.2667 


Spi2= " = — -0. 5334 
1 
q (eqn. 6.7) 


Thus, instead of a matrix of dimension 2” x 2n, we need only 
matrices of dimension (n+l) x (ntl). 

Using MATLAB program we get our new eigenvalues and new 
€igenvectors from aggregated transition matrix Pa. Then we 
get same exact probabilities as given in Appendix A. And 
uSing Equation 5.17 we have approximated probabilities as 


following Equation 6.3. 


O18) 


K ko 
Pp = I - 0.7664 (0.4043) (ean. 6.8) 


0n0 





3p3 


Figure 6.3 Illustration aggregated transition states 


or, Pa: 


We notice that after our first approximation and 
aggregation of the transition matrix P, the probabilities 
are very close. Indeed after certain applications of 
diagnostic Algorithm_1. They are exactly same as seen in 


following tables. 
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TABLE 3 
THE PROBABILITY OF REPAIR OF THE 
(1/2 nodes bad & 2/10 spares bad) SYSTEM 


# OF STEPS 


application of algorithm 1 


1 
1 


POU DONHAUP WN} 


vw 
ve 
vw 
Lh] 
8 
eS 
1 
" 
1 
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TABLE 
THE PROBABILITY OF REPAIR OF THE 
(4/5 nodes bad & 2/10 spares bad) SYSTEM 


EXACT 


APPROX. 


# OF STEPS 


ation of algomithmal 
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VII. CONCLUSIONS AND SUGGESTIONS FOR FUTURE WORK 


A. CONCLUSIONS 

In this thesis, we have considered the effect of faulty 
spares on the probability of repair of digital systems using 
feorithm 1. The analysis utilizes the graph-theoretic 
model of Preparata_Metze_Chien [Ref. 1} and t/s measure of 
Friedman [Ref ©]. 

We realized that for ann node system, we need a 2n x 
on matrix, which is very large for even moderate size n. 
MATLAB can not handle this size of matrix. Then, we get our 
mensey approximation which is Equation 6.3 for two node 
=yovem. 

In Figure 7.1, we have the probability of repair for 
two node system. Indeed, exact and approximate values are 
very close. For example, after first application of 
Algorithm _1, we have an exact value which is 0.6294 and an 
approximate value which is 0.6889. After the second 
application of Algorithm_1, we have an exact value which is 
0.8589 and an approximate value which is 0.8742. Finally, 
after six applications of Algorithm_1, we have the same 
exact and approximate values. 

In Figure 7.2, we have the probability of repair for a 
three node system which has two faulty nodes. After sixteen 


applications of Algorithm_1, the three node system is 
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totally fault-free. The exact and approximate values are 
very close also. After six applications of Algorithm_1; we 
have the same exact and approximate values. 

In Figure 7.3, we have the probability of repair for a 
four node system which has three faulty nodes. In Figure 
7.4, we have the probability of repair for a five node 


system which has four faulty nodes. 


TABLE 7 


# OF NODES IN THE SYSTEM VS 
AGGREGATED TRANSITION MATRIX DIMENSION 


7 OF NODES AGGREGATED TRANSITION MATRIX DIMENSION 
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BEGer that we calculated our aggregated transition 
matrix Pa as a function of the elements of transition matrix 
ie As we see in Table 7 after aggregation, we need only 
matrices of dimension (ntl) x (ntl) instead of a matrix of 
dimension 2” x @2n, 

We calculated the probability of repair for two, three, 
four and five node systems (approximate vs. exact) in 
Meares {/.1, f.2, 7.3 and 7.4. After that we notice our 


first approximation and aggregation of transition matrix P 


the probabilities are almost the same. In fact after 
certain steps of diagnostic Algorithm_1, they are exactly 
the same. 


ae SUGGESTIONS FOR FUTURE WORK 

In this work, we defined an approach which calculates 
each transition state probability. In this thesis, we use 
Algorithm_1 to get probability of repair for two, three, 
four and five node system. We have very close agreements 
with exact solutions. 

One fruitful area of research might be to extend this 
analysis to Algorithm_2 and 3, deriving approximation 
equations for the probability of repair. 

So far we have done many calculations for two node 
systems. Three, four and five node systems need many more 
calculations. Especially, a computer program can be written 
fomtind the transition matrix P. it is hoped that this work 


shows the way for future research. 
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APPENDIX A 
1/2 NODES BAD AND 2/10 SPARES BAD SYSTEM CASE RESULTS 


This eppendix contains matrix manipulations for two node 
system. (Faulty node=1/Faulty spares=2/Fault-free spares=8) 


aE 0.6294 0.6294 0.1666 

0 Ome t2 O. 13a OR 2667 

P= 0 Om eee On geme O 20167 
0 0.0082 0.0082 0.3000 ee 
x 


4 EIGENVECTORS OF THE TRANSITION MATRIX 


a. -0.0000 1.0000 1.0000 
VO=/] O V1i=! 1.0000 V2=|-0.4636 V3=/!-0.6216 
0 -1.0000 -0.4636 -0.6216 
0 -0.0000 -0.0729 0.2431 
4x1 4x1 4x1 4x1 
4 EIGENVALUES OF EO=1.0000 
THE TRANSITION E1=0.1000 
MATRIX E2=0.4043 
E3=0.2581 
INVERSE MATRIX OF EIGENVECTORS 
1 1.0000 1.0000 1.0000 
0 0.5000 -0.5000 -0.0000 
0 -0.7694 -0.7694 —-3.9338 
O -0.2306 -0.2306 2.9338 i 
4x 


TS 


lst Power of Transition matrix P 


2nd Power of Transition matrix P 


OO 
WWW © 
WWW O 
ANNOY 


OOOO 


VAAN 
O40 
NOOO 
WOANO 


OOOO 


VAAN 
O40 
NOOO 
ONAO 


elelele) 


AOOO 


3rd Power of Transition matrix P 


~OOO 
elelele) 


4th Power of Transition matrix 


iE 


SAAN 


OOOO 
OOOO 


| AOOO 


5th Power of Transition matrix 


2 


yoow 


HOOOO 
OOOO 


AOOO 
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6th Power of Transition matrix P 


7th Power of Transition matrix P 


8th Power of Transition matrix P 


9th Power of Transition matrix P 


10th Power of Transition matrix P 
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11th Power of Transition matrix P 


1 1.0000 1.0000 0.9998 | 

0 0.0000 Q.0000 0.0001 

0 0.0000 Q0.0000 0.0001 

0 0.0000 Q0.0000 Q0.0000 Pers 
x 


Al 1.0000 1.0000 Os9999 

0 0.0000 0.0000 Q0.0000 

'O 0.0000 0.0000 Q.0000 

la 0.0000 0.0000 0.0000 on 
x 


13th Power of Transition matrix P 


OO0OF 
OO0OF 


Note : The (1/2 nodes bad & 2/10 spares bad) system is repaired 
after 13 applications of the Algorithm 1. 


ae 


AGGREGATED TRANSITION MATRIX OF THE 2 NODES SYSTEM 


1 0.6294 0.1666 
Pa = 0 0.3624 0.5334 
0 0.0082 0.3000 


3 EIGENVECTORS OF THE AGGREGATED TRANSITION MATRIX 


Ae 1.0000 -0.8044 
VO=| O V1l=|-0.9271 V2=| 1.0000 
0 -0.0729 =O- o> 6 
3x 3x1 3x 
3 EIGENVALUES OF EO=1.0000 
THE AGGREGATED E1=0.4043 
TRANSITION MATRIX BE2=0.2581 


INVERSE MATRIX OF EIGENVECTORS 


OOF 
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lst Power of the aggregated transition matrix Pa 


i 0.6294 0 
0 One 624 0 
0 OSO0GZ 0 


W OIF 
OWd?) 
OWO? 
Of 0) 


2nd Power of the aggregated transition matrix Pa 


~ooG2 oe 

0 Qe ie OE 5533 

0 0054 0944 
3X3 


OO 
Oo 
OoOWw 
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6th Power of the aggregated transition matrix Pa 


8th Power of the aggregated transition matrix Pa 


OOO 


«9994 0.997 2 

-0005 0.0025 
0000 0.0002 

SS 


OOF 


9th Power of the aggregated transition matrix Pa 
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llth Power of the aggregated transition matrix Pa 


ly 1.0000 oo ao 
O 0.0000 0.0001 
0 - 0000 -0000 
3X3 


13th Power of the aggregated transition matrix Pa 


© 
OOF 
OOO 
elete) 
OOO 
OOO 


1.0000 
0.0000 
0.0000 

3X3 


Note : After aggregation the (1/2 nodes bad & 2/10 spares 
bad) system is also repaired after 13 applications 
of the Algorithm 1. 
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APPENDIX B 


2/3 NODES BAD & 2/10 SPARES BAD SYSTEM CASE RESULTS 


This appendix contains matrix manipulations for three node 


system. 


OOOO0O00rF 


(Faulty nodes=2/Faulty spares=2/Fault-free spares=8) 


< 
oO 


OOO0O000rF 





0.5091 
O 


0.3091 
O. foul 


0 
0.0227 
0 
0 


0.5091 0.1712 
0.1591 0.1220 
O 0.4368 
0.3091 0.0565 
QO 0Q.1277 
0 0.0736 
0.0227 Q.0081 
0 0.004) 


8 EIGENVECTORS OF THE TRANSITION MATRIX 





-0.0008 


0.0000 
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0.1712 0.1712 0.0548 
0.0565 0.4368 0.1068 
0.1220 0.0565 0.1068 
0.4368 0.1220 0.1068 
0.0081 0.0736 0.1541 
0.1277 0.0081 0.1541 
0.0736 90.1277 0.1541 
9.0041 0.0041 0.1625 
V3 

0.0000 + 

=O. 50005— 
1.0000 + 

=0- 5000" + 
020153. — 
0.0916 + 

-0.0000 + 
V6 
re 00001 0000 
= 86601 =O 7013 
= 00001 =0 7 ES 
T 86601 =0.7 2 
+ 0.40731 0.4218 
= Soo. 0.4218 
+ 12431 0.4218 
ae 00001 1316 





8x8 


8 EIGENVALUES OF EO=1.0000 


THE TRANSITION E1=0.5150 

MATRIX E2=0.2415 + 0.17681 
E3=0.2415 - 0.17681 
E4=0.2020 
E5=0.0749 + 0.01771 
E6=0.0749 - 0.01771 
E7=0.1231 


INVERSE MATRIX OF EIGENVECTORS 





1 1.0000 - 0.0000i 1.0000 + 0.0000i 1.0000 + 0.0000i 

O -0.9287 + 0.0000i -0.9287 - 0.0000i -0.9287 - 0.00001 

O -0.1659 - 0.22961 0.2818 - 0.0288i -0.1159 + 0.25851 

O -0.1659 + 0.22961 0.2818 + 0.02881 -0.1159 - 0.2585i 

O -0.0484 + 0.0000i -0.0484 - 0.0000i -0.0484 - 0.00001 

O -0.0008 - 0.0590i 0.0515 + 0.0288i -0.0507 + 0.0302i 

O -0.0008 + 0.0590i 0.0515 - 0.02881 -0.0507 - 0.0302i 

O -0.0229 + 0.0000i -0.0229 - 0.0000i -0.0229 - 0.0000i 

1.0000 - 0.0000i 1.0000 + 0.0000i 1.0000 - 0.0000i 1.0000 
-1.9150 + 0.0000i -1.9150 - 0.00001 -1.9150 + 0.00001 -3.3556 

0.3070 — 0.40651 0.1985 + 0.46911 -0.5055 - 0.06271 -0.0000 

0.3070 + 0.40651 0.1985 - 0.46911 -0.5055 + 0.06271 -0.0000 

0.5674 - 0.00001 0.5674 - 0.0000i 0.5674 + 0.00001 6.2433 
-0.3070 + 0.40651 -0.1985 - 0.46911 0.5055 + 0.06271 -0.0000 
-0.3070 - 0.40651 -0.1985 + 0.46911 0.5055 - 0.06271 -0.0000 

0.3475 - 0.00001 0.3475 + 0.0000i 0.3475 + 0.00001 -3.8878 


lst Power of Transition matrix P 
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2nd Power of Transition matrix P 
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3rd Power of Transition matrix P 


OM MH OwWomM 


OAAAOOOO 


OOO9O0000 


MO O10 OM ALN 
MANAOOUWMO 
SPHAOMAOOO 
~OO0O0000 


OOO0O0000 


M010 10 [> tH OLN 


~OO0O0000 


OOOO0000 


MOO HOM LO 
MONDMNOWOS 
TO ADOHOO 
~OOO0000 


OO9O00000 


RdAMNNOdwOd 


NAOOMAMO 
mUNMOOOO 
DOOO90O000 


OOOO0000 


~KOaNATDONe 
NOOSAMMO 
KNMMNOOOO 
DOOO0O000 


OOOO0O000 


~KNAMNOAd 
NO AHOMMAO 
~MNONOOOO 
OODO90O000 


OVOOO0000 


8x8 





| AODOOO000 


4th Power of Transition matrix 
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e 


5th Power of Transition matrix 
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6th Power of Transition matrix P 
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7th Power of Transition matrix P 
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8th Power of Transition matrix P 
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9th Power of Transition matrix P 
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10th Power of Transition matrix P 


Otsdt dAdo 


ANOOOOO0°O 


ODOOO0O000 


NOOO dAdHHAO 


HNOOOO0O00O 


OOOO O000 


MODddaddaAdO 


HNOODOO0O00O 


OOOO0O000 


NOOO AdAO 


OOOOOO00O 


OOOO O000 


ors t OOOO 
MOO OODO000 
OOOOVO0O00 
HNOOO9O0O00O 


OOOOVO0O00 


odds OO0O0 
OCOOOO00°0 
HNOOGVOO90O0 
NOOOOO0O 


OOOO 9O0O00 


ost OOOO 
COODDO0O00O 
HOOOWOOO0O 
HNOOOOO00O 


OOOO O000 


8x8 





| AOOOO0O00 | 


llth Power of Transition matrix P 
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12th Power of Transition matrix P 
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13th Power of Transition matrix P 
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14th Power of Transition matrix P 
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15th Power of Transition matrix P 
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16th Power of Transition matrix P 
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17th Power of Transition matrix P 
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The (2/3 nodes bad & 2/10 spares bad) system is repaired 


Note 


after 17 applications of the Algorithm 1. 
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AGGREGATED TRANSITION MATRIX OF THE 3 NODES SYSTEM 


4 EIGENVECTORS OF THE AGGREGATED TRANSITION MATRIX 


1.0000 
0.5150 
0.2020 
Ooo. 


EO= 


4 EIGENVALUES OF 
THE AGGREGATED 


E1l= 


E2= 


TRANSITION MATRIX 
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INVERSE MATRIX OF EIGENVECTORS 


on) 


Ist Power of the aggregated transition matrix Pa 


2nd Power of the aggregated transition matrix Pa 


3rd Power of the aggregated transition matrix Pa 
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4th Power of the aggregated transition matrix Pa 


5th Power of the aggregated transition matrix Pa 
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6th Power of the aggregated transition matrix Pa 


7th Power of the aggregated transition matrix Pa 


8th Power of the aggregated transition matrix Pa 


9th Power of the aggregated transition matrix Pa 


10th Power of the aggregated transition matrix Pa 


oS 


11th Power of the aggregated transition matrix Pa 


12th Power of the aggregated transition matrix Pa 
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13th Power of the aggregated transition matrix Pa 


14th Power of the aggregated transition matrix Pa 


15th Power of the aggregated transition matrix Pa 
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16th Power of the aggregated transition matrix Pa 


OO0OF 
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4x4 


Note : After aggregation the (2/3 nodes bad & 2/10 spares 
bad) system is also repaired after 17 applications 
of the Algorithm 1. 
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APPENDIX C 
3/4 NODES BAD & 2/10 SPARES BAD SYSTEM CASE RESULTS 
This appendix contains matrix manipulations for four node 
system. (Faulty nodes=3/Faulty spares=2/Fault-free spares=8) 
States 0 through 6 


1.0000 0.4087 0.4087 0.4087 0.4087 0.1099 
O 0.3739 0 0 0.1739 0.3888 
0 O24739 Oes7 359 0 0 0.0547 
0 0 Oc1739 QO. 3750 
0 0 0 0.1739 0.3739 0.1342 
O 0 0 0.0435 0.1842 
0 0 0.0435 0 0 
O 0 0 0 0 0.0182 
P= O 0 0 0 O 
0 0.0435 0 0 0 0.0978 
O 0 0 0.0435 0 O 
0 0 0 0 0 0 
0 0 0 0 O O 
0 0 0 0 O 0.0122 
0 0 0 0 0 0 
0 0 0 0 0 0 
States 7 through 13 
021505 021505 0.1099 0.1099 0.0410 0.0410 
0.0852 0.2194 0.1342 0.0547 0.0204 0.1342 
0.2194 O20352 0.3888 0.0486 0.0204 
0.0852 0.2194 0.0547 0.1342 0.1342 0.0486 
0.2194 0.0852 0.3888 0.1342 0.1342 
0.0365 0.0365 0.0978 0.0553 0.2885 
020365 0.0365 0.0978 0.0702 0.0051 
0.0662 0.0183 0.0702 0.0333 
0.0183 0.0662 0.0182 0.0182 0.0333 0.0702 
0.0365 0.0365 0.1842 0.0051 920322 
0.0365 0.0365 0 0.1842 0.2885 0.0702 
0.0029 0.0020 0 0.0699 0.0041 
0.0020 0.0029 0 0.0122 0.0409 0.0699 
0.0029 0.0020 0 0 0.0041 0.0409 
0.0020 0.0029 0.0122 0 0.0041 0.0041 
0 0.0020 0.0020 


~/ 
ct 
ty 
ct 
(D 
” 
H 
ip 
© 
.Q 
2 F 
RR 
O1 


0.0410 0.0171 
0.0486 0.03387 
0.1342 0.0387 
0.1342 0.0387 
0.0204 0.0387 
070051 0.0651 
0.2885 0.0651 
0.0333 0.0651 
0.0702 0.0651 
0.0702 0.0651 
0.0333 0.0651 
0.0409 0.0875 
0.0041 0.0875 
0.0041 0.0875 
0.0699 0.0875 
0.0020 0.0875 
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lst Power of Transition matrix P 


States O through 6 
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2nd Power of Transition matrix P 


States O through 6 


1.0000 0.6374 0.6374 0.6374 
0 0.1456 0.0326 
0 0.1470 0.1456 0 
0 0.0326 0.1470 0.1456 
0 0 0.0326 0.1470 
0 O 0 0.0118 
0 0.0118 0.0243 0 
0 0 0.0008 0 
0 0.0008 0 0.0008 
0 0.0243 0 0 
0 0.0118 0.0243 
0 0 0.0005 0 
0 O 0 0.0005 
0 O 0 0 
0 0.0005 0 0 
0 0 0 0 
States 7 through 13 
0.4286 0.4286 0.3802 0.3802 
0.1016 0.1349 0.0795 0.1418 
0.1349 0.1016 0.2566 0.0167 
0.1016 0.1349 0.1418 0.0795 
0.1349 0.1016 0.0167 0.2566 
0.0238 0.0181 0.0007 0.0571 
0.0238 0.0181 0.0571 0.0007 
0.0066 0.0042 0.0025 0.0025 
0.0042 0.0066 0.0054 0.0054 
0.0181 0.0238 0.0413 0.0130 
0.0181 0.0238 0.0130 0.0413 
0.0010 0.0009 0.0017 0.0001 
0.0009 0.0010 0.0001 0.0032 
0.0010 0.0009 0.0001 0.0017 
0.0009 0.0010 0.0032 0.0001 
0.0000 0.0000 0.0000 0.0000 


0.2430 0.1444 
0.0586 0.1115 
0.1503 Q0.1115 
0.2275 0.1115 
0.0814 0.1115 
0.0122 0.0652 
0.0929 O. 

0.0146 0.0317 
0.0140 QO. 

0.0259 0.0652 
0.0585 0.0652 
0.0097 0.0192 
0.0031 0.0192 
0.0014 0.0192 
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3rd Power of Transition matrix P 


States O through 6 


1.0000 0.7744 0.7744 0 
0 0.0579 0.0064 0 
O 0.0914 0.0579 0 
0 0.0439 0.0914 O 
0 0.0064 0.0439 0 
0 0.0000 0.0026 0 
0 QO.0111 0.0109 O 
0 0.0002 0.0005 O 
0 0.0005 0.0002 0.0005 
0 0.010 0.0000 0.0026 
0 0.0026 O20liT 0.0109 
0 0.0002 0.0003 0.0000 
0 0.0000 0.0002 0.0003 
0 Q.0000 0.0000 0.0002 
0 0.0003 0.0000 O. 
) 0.0000 0.0000 0. 
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States 7 through 13 
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4th Power of Transition matrix P 


States O through 6 


1.0000 0.8588 0.8588 0.8588 
O 0.0245 0.0118 020578 
O 0.0501 0.0245 0.0118 
0 0.0378 0.0501 0.0245 
0 0.0118 0.0378 0.0501 
0 0.0006 0.0036 O-eo72 
0 0.0072 0.0046 0.0006 
0 0.0003 0.0003 0.0003 
0 0.0003 0.0003 0.0003 
0 0.0046 0.0006 0.0036 
0 0.0036 0.0072 0.0046 
0 0.0002 0.0002 0.0000 
0 ° 9000 0.0002 0.0002 
0 0.0000 0.0000 0.0002 
0 0.0002 0.0000 0.0000 
0 0.0000 0.0000 0.0000 


0.7675 Oe7675 0.7427 0.7427 
0.0502 0.0510 0.0263 0.0845 
0.0510 0.0502 0.0666 0.0454 
0.0502 0.0510 0.0845 0.0263 
0.0510 0.0502 0.0454 0.0666 
0.0072 0.0068 0.0039 0.0124 
0.0072 0.0068 0.0124 0.0039 
0.0007 0.0006 0.0008 0.0008 
0.0006 0.0007 0.0007 0.0007 
0.0068 0.0072 0.0044 0.0113 
0.0068 0.0072 0.0113 0.0044 
0.0002 0.0002 0.0004 0.0001 
0.0002 0.0002 0.0003 0.0002 
0.0002 0.0002 0.0001 0.000 

0.0002 0.0002 0.0002 0.0003 
0.0000 0.0000 0.0000 0.0000 
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5th Power of Transition matrix P 


States O through 6 
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ODDDOC OOOO 0O00000 


0.8535 
-0320 
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ODD OOO0O0O0000000 
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0.9115 0.9115 
0.0124 0.0130 
0.0259 0.0124 
0.0265 0.0259 
0.0130 0.0265 
0.0010 0.0031 
0.0040 0.0020 
0.0002 0.0002 
0.0002 0.9002 
0.0020 0.0010 
0.0031 0.0040 
0.0001 0.0001 
0.0001 0.0001 
0.0000 0.0001 
0.0001 0.0000 
0.0000 0.0000 
0.8376 

0.0207 

0.0333 

0.0501 

0.0376 

0.0040 

0.0058 

0.0004 

0.0004 

0.0024 

-0072 

0.0002 

0.0001 

0.0001 

0.0000 

15 
16x16 
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6th Power of Transition matrix P 
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States 7 through 13 
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7th Power of Transition matrix P 


States O through 6 


1.0000 0.9651 0.9651 O 

0.0056 0.0082 O 

0 0.0071 0.0056 O 

O OmO0CI7 0.0071 0 

O 0.0082 0.0097 0 

O 0.0009 0.0014 O 

O 0.0011 0.0007 O 

0 0.0001 0.0001 0 

0 0.0001 0.0001 O 

O 0.0007 0.0009 O 

O 0.0014 0.0011 O 

0 0.0000 0.0000 O 

O 0.0000 0.0000 0 

O 0.0000 0.0000 O 

O 0.0000 0.0000 0 

0 0.0000 0.0000 0 

States 7 through 13 

0.9421 0.9421 OR 9357 O 

020127 0.0127 0.0127 0 

0.0127 0.0127 0.0110 O 

0.0127 OFO1L27 0.0155 0 

@.0127 0.0127 0.0271 0 

0.0017 0.0017 0.0022 0 

0.0017 0.0017 0.0015 0 

0.0001 0.0001 0.0001 O 

0.0001 0.0001 0.0001 O 

0.0017 0.0017 0.0014 0 

0.0017 0.0017 020023 O 

0.0000 0.0000 0.0000 O 

0.0000 0.0000 0.0001 O 

0.0000 0.0000 0.0000 0 

0.0000 0.0000 0.0000 0 

0.0000 0.0000 0.0000 O 

States 14 through 15 
0.9146 0.8901 
0.0202 0.0240 
0.0159 0.0240 
OF OL/2 0.0240 
0.0215 0.0240 
0.0030 0.0033 
0.0020 0.0033 
0.0002 0.0003 
0.0002 0.0003 
0.0025 020033 
0.0025 0.0033 
0.0000 0.0001 
0.0001 0.0001 
0.0001 0.0001 
0.0001 0.0001 
0.0000 0.0000 
16x16 
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8th Power of Transition matrix P 


States O through 6 


1.0000 0.9781 0.9781 0.9781 
0.0041 0.0055 0.0056 
0 ().0041 0.0041 O20d>5 
0 0.0056 0.0041 0.0041 
0 O.0055 0.0056 0.0041 
0 0.0007 0.0008 0.0006 
0 0.0006 0.0005 0.0007 
0 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 
0 0.0005 0.0007 0.0008 
0 0.0008 0.0006 0.0005 
0 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.9000 
0 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 
States 7 through 13 
0.9636 0.9636 0.9596 0.9596 
0.0080 0.0080 0.0090 0.0087 
0.0080 0.0080 0.0073 0.0104 
0.0080 0.0080 0.0087 0.0090 
0.0080 0.0080 0.0104 0.0073 
QO.0011 0.0011 0.0014 0.0009 
0.0011 0.0011 0.0009 0.0014 
0.0001 Q.0001 0.0001 0.0001 
0.0001 0.0001 0.0001 0.0001 
0.0011 0.0011 0.0011 0.0013 
0.0011 0.0011 0.0013 0.0011 
0.0000 0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 0.0000 
States 14 through 15 
0.9463 0.9308 
0.0130 0.0151 
0.0109 0.0151 
0.0105 0.0151 
0.0126 0 
0.0018 O 
0.0013 O 
0.0001 O 
0.0001 0 
0.0017 0 
0.0014 0 
0.0000 O 
0.0000 0 
0.0000 O 
0.0000 0.0001 
0.0000 0 
16x16 
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9th Power of Transition matrix P 


States O through 6 


1.0000 0.9862 0.9862 O 
0 0.0029 0.0035 O 
0 0.0026 0.0029 O 
0 0.0032 0.0026 0 
0 0.0035 0.0032 0 
0 0.0005 0.0005 0 
O 0.0003 0.0003 0 
0 0.0000 0.0000 QO. 
0 0.0000 0.0000 0 
0 0.0003 0.0005 0 
O 0.0005 0.0003 0 
O 0.0000 0.0000 0 
O 0.0000 0.0000 0 
0 0.0000 0.0000 0 
0 0.0000 0.0000 0 
0 0.0000 0.0000 0 


@.9772 O39772 0.9746 0 
0.0050 0.0050 0.0060 O 
0.0050 0.0050 0.0049 O 
0.0050 0.0050 0.0052 O 
0.0050 0.0050 0.0062 O 
0.0007 0.0007 0.0009 O 
0.0007 0.0007 0.0006 O 
0.0001 0.0001 0.0001 O 
0.0001 0.0001 0.0001 O 
0.0007 0.0007 0.0007 0 
0.0007 0.0007 0.0007 O 
0.0000 0.0000 0.0000 0 
0.0000 0.0000 0.0000 O 
0.0000 0.0000 0.0000 O 
0.0000 0.0000 0.0000 0 
0.0000 0.0000 0.0000 0 
States 14 through 15 

0.9663 O 

0.0081 O. 

0.0073 O 

0.0067 0 

0.0075 0 

0.0010 O 

0.0009 0.0 

0.0001 0 

0.0001 0 

0.0011 0 

0.0009 0 

0.0000 O 

0.0000 O 

0.0000 0 

0.0000 O 

0.0000 0 

16x16 
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OOO OOO0O0O00000000 
OOO DOO0O000000000 


OOOO ODOCO OO 0000000 
OOVDODOOOOO 0000000 


10th Power of Transition matrix P 


States O through 6 


1.0000 0.9914 0.9914 0.9914 
O 0.0019 0.0021 0.0019 
0 0.0017 0.0019 0.0021 
0 0.0019 0.0017 0.0019 
0 0.0021 0.0019 0.0017 
0 0.0003 0.0003 0.0002 
O 0.0002 0.0002 0.0003 
0 0.0000 0.0000 0.0000 
O 0.0000 0.0000 0.0000 
O 0.0002 0.0003 0.0003 
0 0.0003 0.0002 0.0002 
O 0.0000 0.0000 0.0000 
O 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 
O 0.0000 0.0000 0.0000 

States 7 through 13 

0.9857 0.9841 0.9841 

0.0031 0.0038 0.0032 

0.0031 0.0033 0.0037 

0.0031 0.0032 0.0038 

0.0031 0.0037 0.0033 

0.0004 0.0005 0.0004 

0.0004 0.0004 0.0005 

0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 

0.0004 0.0005 0.0004 

0.0004 0.0004 0.0005 

0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 

States 14 eugn 15 

0.9788 

0.0050 

0.0048 

0.0043 

0.0045 

0.0006 

0.0006 

0.0000 

0.0000 

0.0007 

0.0006 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


16x16 
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llth Power of Transition matrix P 
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12th Power of Transition matrix P 


States O through 6 


1.0000 G229266 0.9966 O 
O 0.0008 0.0008 O 
O 0.0007 O20003 0 
O 0.0007 0.0007 O 
0 0.0008 0.0007 O 
O 0.0001 0.0001 O 
O 0.0001 0.0001 O 
O 0.0000 0.0000 O 
O 0.0000 0.0000 0.0000 
O 0.0001 0.0001 O 
O 0.0001 0.0001 0 
0 0.0000 0.0000 O. 
0 0.0000 0.0000 GO: 
0 0.0000 0.0000 0 
0 0.0000 0.0000 0 
0 0.0000 0.0000 OF 
States 7 through 13 
0.9943 0.9937 O 
0.0012 0.0014 O 
0.0012 0.0014 0 
0.0012 0.0013 0 
0.0012 0.0013 O. 
0.0002 0.0002 0 
0.0002 0.0002 0 
0.0000 0.0000 0.0000 
0.0000 0.0000 O. 
0.0002 0.0002 0 
0.0002 0.0002 0 
0.0000 0.0000 0 
0.0000 0.0000 0 
0.0000 0.0000 Oe 
0.0000 0.0000 OF 
0.0000 0.0000 0.0000 
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13th Power of Transition matrix P 


States O through 6 


1.0000 0.9979 0.9979 QO. 
0 0.0005 0.0005 0.0004 
0 0.0005 0.0005 0.0005 
0 0.0004 0.0005 0.0005 
0 0.0005 0.0004 0 
0 0.0001 0.0001 0 
0 0.0001 0.0001 0 
O 0.0000 0.0000 0 
0 0.0000 0.0000 0.0000 
0 0.0001 0.0001 0 
0 0.0001 0.0001 0 
0 0.0000 0.0000 0 
0 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0. 
0 0.0000 0.0000 QO. 





0.9964 0.9964 0.9961 0.9961 
0.0008 0.0008 0.0009 0.0008 
0.0008 0.0008 0.0009 0.0008 
0.0008 0.0008 0.0008 0.0009 
0.0008 0.0008 0.0008 0.0009 
0.0001 0.0001 0.0001 0.0001 
0.0001 0.0001 0.0001 0.0001 
0.0000 0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 0.0000 
0.0001 0.0001 0.0001 0.0001 
0.0001 0.0001 0.0001 0.0001 
0.0000 0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 Q.0000 
0.0000 0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 O. 
States 14 through 15 

0.9948 OFI9S2 

0.0011 0.0015 

0.0012 0.0015 

0.0012 0.0015 

0.0011 0.0015 

0.0001 0.0002 

0.0002 0.0002 

0.0000 0.0000 

0.0000 0.0000 

0.0002 0.0002 

0.0002 0.0002 

0.0000 0.0000 

0.0000 0.0000 

0.0000 0.0000 

0.0000 0.0000 

0.0000 0.0000 

| 16x16 


109 


14th Power of Transition matrix P 


States O through 6 


OODOO0000000 
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1.0000 0.9987 0.9987 0.9987 
O 0.0003 0.0003 0.0003 

O 0.0003 0.0003 0.0003 

O 0.0003 0.0003 0.0003 

O 0.0003 0.0003 0.0003 

O 0.0000 0.0000 0.0000 

0 O0.0000 0.0000 0.0000 

0 0.0000 0.0000 0.0000 

0 0.0000 0.0000 0.0000 

O 0.0000 0.0000 0.0000 

O 0.0000 0.0000 0.0000 

O 0.0000 0.0000 0.0000 

0 0.0000 0.0000 0.0000 

J 0.0000 0.0000 0.0000 

O 0.0000 0.0000 0.0000 

0 0.0000 0.0000 0.0000 

States 7 through 13 
0.9978 0.9978 O2397> 0.9975 
0.0005 0.0005 0.0005 0.0005 
0.0005 0.0005 0.0006 0.0005 
0.0005 0.0005 0.0005 0.0005 
0.0005 0.0005 0.0005 0.0006 
0.0001 0.0001 O.0001 0.0001 
0.0001 0.0001 0.0001 0.00901 
0.0000 0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 0.0000 
0.0001 0.0001 0.0001 0.0001 
0.0001 0.0001 0.0001 0.0001 
0.0000 0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.0000 
15 
16x16 





110 


15th Power of Transition matrix P 


States O through 6 


1.0000 0.9992 0.9992 On 9292 
O 0.0002 0.0002 0.0002 
O 0.0002 0.0002 0.0002 
O 0.0002 0.0002 0.0002 
O 0.0002 0.0002 0.0002 
O 0.0000 Q0.0000 0.0000 
0 0.0000 0.0000 0.0000 
O 0.0000 0.0000 0.0000 
O 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 
O O.0000 0.0000 0.0000 
0 Q0.0000 QO.0000 0.0000 
O QO.0000 0.0000 0.0000 
O QO.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 


0.9986 0.9984 0.9984 
0.0003 0.0003 0.0003 
0.0003 0.0003 0.0003 
0.0003 0.0003 0.0003 
0.0003 0.0003 0.0003 
0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 
States 14 15 
0.9979 
0.0004 
0.0005 
0.0005 
0.0005 
0.0001 
0.0001 
0.0000 
0.0000 
0.0001 
0.0001 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
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16th Power of Transition matrix P 


States O through 6 


1.0000 Of2295 O23995 0 
0.0001 0.0001 0.0001 
0 0.0001 0.0001 0 
O 0.0001 0.0001 0 
0 0.0001 0.0001 0.0001 
O 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 
O 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 
O 0.0000 0.0000 0.0000 
0 0.0000 0.0000 QO. 
O 0.0000 0.0000 0 
0 0.0000 0.0000 0 
0 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0 
0 0.0000 0.0000 O. 
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17th Power of Transition matrix P 


States O through 6 
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18th Power of Transition matrix P 


States 0 through 6 


1.0000 0.9998 0.9998 0.9998 0.9998 0.9996 0.9996 
O 0.0000 0.0000 0.0000 0.0000 0.0001 0.0001 
O 0.0000 0.0000 0.0000 0.0000 0.0001 0.0001 
O 0.0000 0.0000 0.0000 0.0000 0.0001 0.0001 
0 0.0000 0.0000 0.0000 0.0000 0.0001 0.0001 
0 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.0000 Q.0000 0.0000 0.0000 0.0000 
O 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.00LyU 0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
States 7 through 13 
029997 0.9997 0.9996 0.9996 0.9995 0.9995 0.9995 
0.0001 Q.0001 0.0001 0.0001 0.0001 0.0001 0.0001 
0.0001 0.0001 0.0001 0.0001 0.0001 0.0001 0.0001 
0.0001 0.0001 0.0001 0.0001 0.0001 0.0001 0.0001 
0.0001 0.0001 0.0001 0.0001 0.0001 0.0001 0.0001 
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
0.0000 0.0000 Q.0000 Q.0000 0.0000 Q.0C00 0.0000 
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 0.000 0.0000 0.000 0.000 
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 Q.0000 0.0000 - 0000 0.0000 
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
0.0000 Q.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 


States 14 through 15 





029995 0.9993 
0.0001 0.0001 
0.0001 0.0001 
Q.0001 0.0001 
0.0001 0.0001 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 


16x16 
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19th Power of Transition matrix P 
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States 7 through 13 
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0.9999 Omo3o9 
0.0000 0.0000 
0.0000 0.0000 
Q0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.9998 0.9998 
0.0000 0.0001 
O0.0000 0.0001 
0.0000 0.0001 
0.0000 0.0001 
O.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
through 15 

0 

0 

0 

O 

0 

O 

0.0000 

0.0000 

0.0000 

0.0000 

0.000 

O 

O 

O 

0 
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20th Power of Transition matrix P 


States O through 6 


1.0000 OF 2999 G9 292 0 
0 0.0000 0.0000 ) 
0 0.0000 0.0000 O 
0 0.0000 0.0000 0 
0 0.0000 0.0000 O 
0 0.0000 0.0000 0 
0 0.0000 0.0000 0 
0 0.0000 0.0000 0.0000 
O 0.0000 0.0000 0 
0 0.0000 0.0000 0 
0 0.0000 0.0000 QO. 
O 0.0000 0.0000 0 
0 0.0000 0.0000 0.0 
0 0.0000 0.0000 0 
O 0.0000 0.0000 0 
0 0.0000 0.0000 0 


States 7 through 13 


0.9999 022999 0239935 0 
0.0000 0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 O 
0.Q000 0.0000 0.0000 O 
0.0000 0.0000 0.0000 O. 
0.0000 0.0000 0.0000 0.0000 
0.0000 0.0900 0.0000 QO. 
0.0000 0.0000 0.0000 0 
0.0000 0.0000 0.0000 0 
0.0000 0.0000 0.0000 0 
0.0000 0.0000 0.0000 0 
0.0000 0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 0 
0.0000 0.0000 0.0000 0 
0.0000 0.0000 0.0000 0 
0.0000 0.0000 0.0000 QO. 
States 14 through 15 

0.9998 Oro 27 

0.0000 0.0001 

0.0000 0.0001 

0.0000 ©) 

0.0000 0 

0.0000 O 

0.0000 0 

0.0000 0 

0.0000 0 

0.0000 O 

0.0000 O 

0.0000 O 

0.0000 0 

0.0000 0 

0.0000 O 

0.0000 0 
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ODODDD OOOO O0O00000 


21st Power of Transition matrix P 


States O through 6 


1.9000 1.0000 1.0000 1.0000 1.0000 Oeo299 Oe ogo 
O 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 0.0000 Q.0000 0.0000 
O 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
O 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
O 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
O 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
O 0.000 0.0000 0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
O 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.0000 90.0000 0.0000 - 0000 QO. O 
O 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
O 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 0.0000 Q.0000 0.0000 

States 7 through 13 

0.9999 0.9999 0.9999 0.9999 0.9999 0.9999 0.9999 

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 Q0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 0.0000 Q0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

States 14 through 15 

0.9999 

Q.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 
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22nd Power of Transition matrix P 


States O through 6 


1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 
0 0.0000 0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 0.0000 0.0000 
O 0.0000 0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 0.0000 0.0000 
0 0.0000 0.0000 0.0000 0.0000 0.0000 

States 7 through 13 

1.0000 1.0000 1.0000 1.0000 0.9999 022925 

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 OF 0 QO. 0 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 

States 14 through 15 

0.9999 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 
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23rd Power of Transition matrix P 


States O through 6 
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OO 
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i 
0.0000 0.0000 0.0000 0 
0.0000 0.0000 0.0000 0 
0.0000 0.0000 0.0000 O 
0.0000 0.0000 0.0000 0 
0.0000 0.0000 0.0000 0 
0.0000 0.0000 0.0000 0 
0.0000 0.0000 0.0000 0 
0.0000 0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 0 
0.0000 0.0000 0.0000 0 
0.0000 0.0000 0.0000 0 
0.0000 0.0000 0.0000 0 
0.0000 0.0000 0.0000 0 
0.0000 0.0000 0.0000 O 
0.0000 0.0000 0.0000 0. 
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1.0000 1.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
0.0000 0.0000 
1.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
ough 15 
16x16 
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Note : The (3/4 nodes bad & 2/10 spares bad) system is repaired 
after 24 applications of the Algorithm 1. 
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1st Power of the aggregated transition matrix Pa 


2nd Power of the aggregated transition matrix Pa 
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llth Power of the aggregated transition matrix Pa 
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13th Power of the aggregated transition matrix Pa 
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16th Power of the aggregated transition matrix Pa 
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22nd Power of the aggregated transition matrix Pa 
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After aggregation the (3/4 nodes bad & 2/10 spares 
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Note 


bad) system is also repaired after 24 applications 


of the Algorithm 1. 
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APPENDIX D 
4/5 NODES BAD & 2/10 SPARES BAD SYSTEM CASE RESULTS 


This appendix contains matrix manipulations for five node 
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lst Power of the aggregated transition matrix Pa 
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26th Power of the aggregated transition matrix Pa 
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31st Power of the aggregated transition matrix Pa 
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Note 


bad) system is also repaired after 34 applications 


of the Algorithm 1. 
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